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ABSTRACT 

A pole placement edgorithm is proposed vhich uses constrained non-linear 
optimization techniques on a finite dimensioned model of a linear n degree of 
freedom system. Low order feedback control is assumed where r poles may be 
assigned; r being the rank of the sensor coefficient matrix. It is shown that by 
combining feedback control theory methods with optimization techniques, one can 
ensure the stability characteristics of a system, emd can alter its trrinsient 
response. 


INTRODUCTION 

One common method of approaching the problems of controlling the vibration 
of a structure is o employ eigenvalue (pole) placement methods. Such solutions 
have attracted the attention of numerous authors over the past twenty-five years, 
including W. M. Wonham [6], E. J. Davison [3], S. Srinathkumar [5l» A. 5. Andry 
et al [l], [2] and many others. 

In exploring pole placement in dynamical systems, em inadequacy of stability 
considerations in contemporary algorithms was noted and thus motivated this work. 
It appears that the problem has not been solved or even addressed in many 
approaches . 

If a system is controllable, one has the ability to place a predetermined 
number of poles. Thus, when pole placement techniques are employed, there is a 
limit on the number of poles that may be assigned. As is well known, the rank of 
the sensor coefficient matrix determines how many poles may be placed exactly. 
These poles may be noted as the contrcllaDle eigenvalues of the system, while the 
remaining may be labelled uncontrollable. 
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Thus, due to restrictions inherent to every system, every pole may not be 
desirably placed. Therefore, one does not have control over the full order of 
the system. When moving the allowable eigenvalues, those which axe not placed 
will also be affected, with the possibility of generating an unstable state. 

Since an unstable system is undesirable, the ability to place a pre- 
determined nvimber of poles, while forcing the system to remain stable woxild be 
quite desirable to the designer. Many pole placement methods yield satisfac- 
tory assignment of the desired modes, but unfortunately caji drive the remaining 
eigenvalues unstable. Thus, requiring iteration of the algorithms, compromising 
the desired choice of eigenvalues or eigenvectors, until a stable response 
results. With the large number of modes required in modelling flexible struc- 
tures, these methods become costly €ind time consuming. 

Hence, a pole placement method is proposed which constrains the unspecified 
modes to be stable by taking advantage of constrained qptimization techniques. 

It appears that no previous work has guaranteed stable unplaced poles or has 
assured the magnitude of relative stability. 

Several numerical examples will be presented, and results will be compared 
with those of Srinathkumar [ 5 ] • 


PROPOSED SOLUTION 

Hie systems studied in this paper are of the mechanical type, which are 
second order by nature, incorporating mass, stiffness tind damping parameters, 
where only the class of discrete systems shall be investigated. 

Assuming small motions about the equillibrium point inq>lies linearization 
of the equations of motion, which become 

[M]q(t) + [D+Gl4(t) + [S+H]£(t) = F(t) (l) 

The forcing function vector, F(t), may then be described as 

F(t) = [V]q(t) + [P]q(t), 

where [v] and [P] are the velocity and position feedback matrices, respectively. 
q(t) is the coordinate •vector, •while ^(t) and q(t) axe the first euid second time 
derivatives of this vector. 

[M] is known as the mass or inertia matrix, [D] is called the damping matrix, 
and [s] is the stiffness matrix. The matrix [G] may be referred to as the gyro- 
scopic or Coriolos matrix, and [H] is the circulatory matrix. 

The [M] , [D], [S], [g] and [H] matrices are assumed to be time-invariant, 
and therefore are represented by constant •values, all be’ng of nth order, where 
n represents the niunber of degrees of freedom of the system. 

Using normal state space methods by letting 
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More simply, equation (2) may te expressed as follows: 
x(t) = [A']x^(t) + [B]u(t), x(0) « 

Z(t) = [C]x(t) 
vhere u(t) = [K]x(t) 

X(t) is the output vector, [CJ is a constant sensor coefficient matrix, and [K] 
is the feedback gain matrix. may now be described as the constant coeffi- 

cient matrix of actuator dynamics, and u(t) is the control vector. The following 
conditions hold: 

i) xeR^*^, ueR°, ^eR** 

ii) A', B, C are rea.1., constemt matrices of appropriate dimensions, 

iii) rank B = m # 0, remk C = r # 0 


By block diagram representation, the system described by equation (3) may be 
expressed as in Figure 1. 



FIGURE I 
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And a more revealing representation is shown in Figure 2. 



FIGURE 2 


Equation (2) may be rewritten as follows: 


x(t) = 

-M"^(D+G) I -M"^(S+H) 

x(t) + 

®1 


I ! 0 

_ n 1 


3. 


[K][C^ I C2lx(t) 


or 


c(t) = 



-M"^(S+H) 

x(t) + 

^^^1 1 2 

I 

0 

B„KC, i B„KC„ 

L n 1 J 


_ 2 11 2 2_ 


x(t) 


£(t) = [c^j CgWt) 

By comparison of equations (l) and (3), one may note that this implies: 

iBg] = [0], 


(3) 


thus 


x(t) = 



-M”^(D+G) ! -M"^(S+H) 

x(t ) + 

Yli 



_ 1 ^ _ 


0 

0 


x(t) 


2;(t) = [C^ jC2]x(t) 


(U) 


If we define 


[A] = 


~M"^(D+G) + B^KC^ 


-M"^(S+H) + B^KC^ 
0 


_l 
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and describe equation (U) as follows: 

£ = [A]x 
Z = [C]x 

Then, the set of equations must satisfy the eigenvalue problem, i.e., 

[*lVi - tjVj 


( 5 ) 


where 


i-i: 


the 2n eigenvalues 


2n 


and 


|v^^^ = the corresponding eigenvectors. 

. By substitution of equation ( 4 ) into equation (5)» 




-M"^(D+G) j -M~^(S+H) 

1 

T ! 0 



B,KC, 

iB,KC« 

V 4 - 

1 1 

1 2 

^ * 

0 

1 


n I 

V. may then be defined to correspond to the above peurtitioning as follows: 


yielding 

% 

which implies 
substituting. 


— — 

z. 

“T. 


w. 

1 _ 



I 


V, = 


h 


-m'^^Cd+g) I -m"^(s+h) 
^n I ^ 


?iVi = 




2i 


B^Ci jB^KC^ 






= -M"^(D+G)C^2£^ - M"^(S+H)Wj + B^KC^w^ 

If we define {X^} i = l,2,..,,r as the r eigenvalues to be placed, 
equation (6) may be expressed as 

WA^ = -M"^(I>4-G)WA - M“^(S+H)W + B^KC^WA + B^KC^W 


( 6 ) 
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where 


I w ] 

I "T 


and A = diag(X ) 

X 4i r 

By taking advan** ^.ge of the generalized left inverse theorem, 

[K] = [B^B^]"^[B^][WA^4M“^(D+G)WA + M“^(S+H)W][C^WA+C2W]"^, 

which is the equation describing the gain matrix needed to obtain those eigen- 
values desired. 

A single objective function was then determined from the set of equations 
described by equt-tion (7)» where the values of [K] were determined by minimising 
that objective function. The constraints imposed on the system were that the 
real part of the eigenvalues of the closed loop system were all negative. These 
constraints were also modified, as was desired, to inc: ease the stability *xargin. 

NUMERICAL E5CAMPLES 


Example 1: 



FIGURE 3 

Specifications r = nig = 1 
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Eigenvalues of unforced system: 


-1.666207 ± l.l+133Ui 
- .333783 ± -.832651 


[C] 


0 ■’ 0 0 
_1 J 0 0_ 


[B] 


1 0 
0 1 
0 0 
0 0 


Desired eigenvalues: 




1 


X 


2 


-4.0 + Oi 
-3.0 + Oi 


Resulting eigenvalues using the proposed method with no additional factor for 
relative stability: 

-4.000000 + Oi 
-1.335420 + Oi 
- .249608 + Oi 

000000 + Oi 


Resulting eigenvalues using the proposed method with added factor of relative 
stability: 

-4.000000 + Oi 
“2.999999 + Oi 

- .572543 ± 0.743531 


Resulting eigenvalues using Srinathkumar method: 

9.1256 + Oi 
- .8l4l + Oi 
-4.0000 + Oi 
-3.0000 + Oi 

Note that the method proposed here yields the vdesired .eigenvalues and that 
the unspecified eigenvalues remain stable, whereas .ca the- SrinatiJtumar method 
an unspecified eigenvalue is moved into the right half plane. 
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Example 2: 



Eigenvalues of unforced system: 


-.OOU055 ± 1. 61+79531 
-.17061+9 ± 1.13ll+l8i 

-.062361+ ± .35567^1 
-.0751+32 ± .7301+l+li 

r 1 _ r 1 0 0 0 0 0 0 0 1 

^ ^ ■ [_0 0 0 1 0 0 0 OJ 


"100 
0 10 
0 0 1 
[B] = 000 

0 0 0 
0 0 0 
_0 0 0 _ 

Desired eigenvalues: o “ - ‘5 
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Resulting eigenvalues using the proposed method, \diere a factor for relative 
stability was added: 


-.28931*2 ± 1.378583i 
-.U5U51 ± 1.1713U5i 
-.1*00007 ± .5000031 
-.19781*0 ± .42591*1*1 

Exeimple 3: 



FIGURE 5 


Specifications: ® = 1 

51 = 3 

52 = l 

Eigenvalues of \inforced system: 

±2.0? 43131 
± .8350001 

[C] =[1100] 

~1 0 ~ 

fs! = 0 0 

0 0 


Desired eigenvalue: .5 + Oi 


Resulting - igenvalues using the propored method, where factor for relative 
stability was added: 


- .170373 ± 1.8090971 
-1.817157 + 01 

- .500000 + 01 


CONCLUSION 

A pole placement algorithm has been proposed which used constrained non- 
linear programming techniques for a finite dimensional model of a linear n degree 
of freedom system. It has been shewn that by constraining the eigenvalues of the 
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fMll order system while simultaneously placing those allowable > one can ensure 
the stability characteristics of a system, and can alter its transient response. 

Results of the Srinathkumar method were presented for Example 1, and showed 
how this method yielded the desired eigenvalues quite accurately, yet unfortu- 
nately forced the originally stable system unstable, therefore resulting in an 
undesirable response. 

No previous work has guaranteed stable unplaced poles or has assured the 
magnitude of relative stability. 
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